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Abstract 

A modification is proposed for the formula known from the literature that characterizes the boundary 
of the capacity region of Gaussian multiaccess fading channels. The modified version takes into account 
potentially negative arguments of the cumulated density function that would affect the accuracy of the 
numerical capacity results. 
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I. Introduction 

The boundary of the capacity region of multiaccess (MAC) fading channels was first 
characterized in [1] and discussed in full detail in [2]. It is assumed that the fading processes of 
all users are independent of each other, are stationary and have continuous probability density 
functions, fi(h) Vi, with h > the random fading coefficient and i the user index; a total of 
M users are assumed. The cumulated density functions of the fading processes are denoted by 
Fi(h) = fi(h')dh'. Note that, according to the standard fading channel model with coherent 
detection, the support of the channel coefficients does not contain negative numbers. The receiver 
noise is assumed to be Gaussian with the variance a 2 . 

II. Boundary of the Capacity Region and Modification of the Standard Result 

It was shown in [2, Theorem 3.16] that the boundary of the capacity region of the Gaussian 
multiaccess channel is the closure of the parametrically defined surface 

|r(//) :petf{,Y ( iH = l^ (1) 

where for each i — M 
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The vector // = {0 < jii < 1 : i = 1,2,..., M} is a given "rate award" vector that is specified 
to pick a desired point on the boundary of the capacity region. The vector A = {A« G 9?+ : i = 
1,2,..., M} is the solution of the equations 

OO s OO >. 
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Mi =x 

(3) 

where Fj is the long-term average power constraint of user i. The solution of (3) for the vector 
A is unique, and an iterative numerical procedure is given in [2] to find it. 

As < < 1 W, the differences fj, k — ^ in (2) and (3) can have negative values and, 
hence, the arguments of the cumulated density functions (CDFs) can, depending on the channel 
coefficient h, also be negative. As the fading coefficients can not be negative, the CDF is actually 
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not defined for such values as they lie outside the support of the random variable. Although it 
seems natural to assume the value "zero" in those cases, which might implicitly happen in a 
implementation of (2) and (3), this would lead to incorrect results as we show below. 

To compensate for this problem, we propose to introduce a modified argument in the cumulated 
density functions F k (x) in the expressions in (2) and (3) as follows: 



F k (x) ^ F k ([x]*) (4) 



with 



and 



x = 2X ^ 2 + *) (5 ) 

2\i(a 2 + z) + (fJ, k ~ fi>i)h 

, x if x > 
[*]*=< ~ • (6) 

+oo if x < 

For negative arguments, x, the function [x\* takes on the value +oo which is inserted into a 
CDF in (4). Hence the value of the CDF for x < is "1" and not "0". The justification is given 
in Section III. 

III. Explanation 

There is no need to go through the whole derivation again to characterize the capacity boundary 
surface. We start at the point where we propose a modification, i.e., equation (18) on page 2804 
of [2]. We wish to compute the rate 



with 



\i, z) = Pr (ui(z) > Uj{z) Vj and Ui (z) > o) (8) 



where the marginal utilities ("rate revenue minus power cost" [2, p. 2802]) are defined by 

To solve (7) (and also the corresponding problem in [2, equation (18)] for the vector A to fulfil 
the average power constraint for the user i) we need to evaluate the probability (8). 
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Firstly, it should be noted that the condition Ui(z) > Uj(z) Vj in (8) (implicitly) excludes the 
case j — % because otherwise P(i, z) would be "zero" as, trivially, P(ui(z) > Ui(z)) = 0. Using 
(9) we can state the equivalence 

Ui (z) > ^ hi > 2X ^ a2 + z ^ > o . (10) 

fJ'i 

Note that Aj > Vi, as A is a Lagrange multiplier that introduces the "power price" (that can 
never be negative) into the optimisation problem that must be solved to find the capacity region 
[2]. 

Using (10), the probability (8) can now be written as 

P(i, z) = Pr (ui(z) > | Ui (z) > Uj (z) Vj) • Pr (u t {z) > Uj (z) Vj) (11) 
= Pr (hi > 2X ^ 2 + Z ^ | Ui ( z ) > u .( z ) Vj ) ■ Pr ( Ui (z) > Uj (z) Vj) (12) 

oo 

= J fi(h\ Ui(z) > Uj (z) Vj)dh • Pr (ui(z) > Uj (z) Vj) (13) 

2^(<t 2 + Z ) 
Mi 

OO 

= J fi(h,Ui(z)> Uj (z)Vj)dh (14) 



Hi 
OO 



= J fi(h) -Pr (ui(z)> Uj (z)Vj\hi = hyh (15) 

H 

Since the fading processes of the users are assumed to be independent, we can write: 

r°° _ 

P(i,z)= f2 f i (h)-l\Pr(u i (z)>u k (z)\hi = h)dh. (16) 

Mi kj±i 

Now, we need to evaluate the probability 

Pr (ui(z) > u k (z) \ hi — h) (17) 

We use (9) to rewrite the event Ui{z) > u k (z) and obtain 

/ \ . / \ , Aj M fe A fc 

Ui{z) > u k (z) — > — (18) 

a hi a h k 

or, equivalently, 

a\ k hi h k 
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with the abbreviation a = 2(<r 2 + z) > and \ > Vz and < < 1 Vi. As /i^ - ^ can 
be negative, the left-hand side of (19) can be negative so we have to differentiate between two 
cases: 



Case A : hi(fj, k — Hi) + K a > 
Case B : hi(fj, k — Hi) + A«a < 



(to > Hi) or {^Hk < Hi and ^ < J^jl^j 



fik < Hi and hi > 



a) Case A: With a = 2(a 2 + 2) we obtain from (17), (19) and (20) 

2\ k h l (a 2 + z) 



Pr (^(z) > u k (z) \h { = h) = Pr yi k < 

= Pr(h k < 

= fJ 



2Aj(a 2 + z) + (/i fe - 
2A fc /t(a 2 + z) 
2\i(a 2 + z) + (fj, k - fii)h 
2X k h(a 2 + z) 



hi = h 



(20) 
(21) 

(22) 
(23) 
(24) 



2\i(a 2 + z) + (/i k - fii)h 
with = Jq f k (h)dh the cumulated density function of the channel coefficient k. The 

solution (24) is the one originally used in equations (2) and (3) that are taken from [2]. 
b) Case B: For a negative left-hand side in (19) we obtain 



with 



Pr ( Ui (z) > u k {z) \hi = h) = Pr (h k > B) = 1 
2X k h(a 2 + z) 



B 



< . 



(25) 



(26) 



2Aj(cr 2 + z) + {Hk - Hi)h 
As h k is a channel coefficient and non-negative by definition, the probability (25) is simply 



c) New formulation of the boundary of the capacity region: In order to keep the structure 
of the original solution given in [2] but with the correct evaluation of the probability in both 
cases A and B, we write the probability 

2X k h(a 2 + z) 



Pr (ui(z) > u k (z) \hi = h) = F k 



(27) 



2Aj(<7 2 + z) + {Hk - Hi)h 
with the function [x\* defined in (6). When we use (27) in (16) and (7) we obtain the corrected 
solution proposed in Section EL 
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